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ABSTRACT 

The masses of central massive black holes in BL Lac objects are estimated from their host galaxy absolute 
magnitude Mr at /?-band by using the empirical relation between absolute magnitude of host galaxy Mr and black 
hole mass M D h- Only a small fraction of BL Lac objects exhibit weak broad-line emission, and we derive the 
sizes of the broad-line regions (BLRs) in these BL Lac objects from the widths of their broad emission lines on 
the assumption of the clouds being virilized in BLRs. It is found that the sizes of the BLRs in these sources are 
usually 2-3 orders of magnitude larger than that expected by the empirical correlation ^blr - AL A (3000A) defined 
by a sample of Seyfert galaxies and quasars (McLure & Jarvis 2002). We discuss a variety of possibilities and 
suggest it may probably be attributed to anisotropic motion of the BLR clouds in these BL Lac objects. If the BLR 
geometry of these sources is disk-like, the viewing angles between the axis and the line of sight are in the range 
of ~ 2° - 12°, which is consistent with the unification schemes. 

Subject headings: galaxies: active — BL Lacertae objects: general — accretion, accretion disks — black hole 
physics 
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1. INTRODUCTION 

There is a tight correlation between the black hole mass Mbh 
and the stellar dispersion velocity a (Ferrarese & Merritt 2000; 
Gebhardt et al. 2000). This tight correlation Mbh - cr is widely 
used to estimate the central black hole masses of active galactic 
nuclei (AGNs). Unfortunately, the stellar dispersion velocity a 
is available only for a small fraction of AGNs. McLure & Dun- 
lop (2002) derived a very tight correlation between host galaxy 
absolute magnitude Mr at /?-band and black hole mass M D h- 
Cao (2003) used this relation Mbh -Mr to estimate the central 
black hole masses of 29 BL Lac objects. The central black hole 
masses of three sources in their sample have also been mea- 
sured from the stellar dispersion velocity (Falomo, Kotilainen, 
& Treves 2002; Barth, Ho, & Sargent 2003), which agree well 
with the black hole masses estimated from the host galaxy abso- 
lute magnitude Mr. O'Dowd, Urry, & Scarpa (2002) also found 
that the black hole masses estimated from the host galaxy lumi- 
nosity are quite reliable for radio galaxies. 

The BLR sizes Rblr of broad-line H/3 are measured by Kaspi 
et al. (2000) for a sample of quasars and Seyfert galaxies using 
reverberation mapping method. They found a tight correlation 
between the BLR size Rblr and optical continuum luminosity 
\L\ . Using the width of the broad emission line and measured 
BLR size Rblr, they estimated the central black hole masses 
of the sources in their sample assuming the clouds in BLRs 
to be virilized. For sources at high redshifts, the emission of 
line H/3 is usually unavailable. Instead, the width of the line 
Mg II can be used to estimate the central black hole masses 
(McLure & Jarvis 2002). Mg II is a low-ionization line as H/3, 
so that Mg II is expected to be produced in the same region as 
H/3, which is supported by the tight correlation between the 
FWHM of Mgll and H/3 (VpwHM(Mg n)~ V F whm(H/3)) found 
by McLure & Jarvis (2002). It is therefore reasonable to ex- 
pect that Mg II and H/3 are produced in the same region. Using 
the same sample and the BLR sizes measured by Kaspi et al. 
(2000), McLure & Jarvis (2002) obtained a correlation between 
Rblr and the monochromatic continuum luminosity at 3000 A, 
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which is useful for estimate of black hole masses of the sources 
at high redshifts with only Mg II emission line profiles. 

For most AGNs, their BLR sizes have not been measured 
directly by reverberation mapping method and the empirical re- 
lation Rblr - ALa is used to derive Rblr- Combining the line 
width, the central black hole masses can be estimated by assum- 
ing the motion of clouds in BLRs to be virilized. The estimated 
black hole masses depend sensitively on the velocity of the 
clouds in BLRs (oc Vb LR ). The broad-line width is mainly gov- 
erned by the component of the cloud velocity Vblr projected 
to the line of sight. If the motion of BLR clouds is anisotropic 
(e.g., disk-like BLR geometry), the estimate of black hole mass 
becomes complicated (e.g., Jarvis & McLure, 2002). McLure 
& Dunlop (2001) argued that BLRs in some AGNs have disk- 
like geometry. If the disk-like BLR geometry is indeed present, 
the observed broad-line width depends sensitively on the orien- 
tation of disk axis and the orientation o is therefore crucial in 
the estimate of black hole mass from its broad-line width. In the 
unification schemes, the jets of BL Lac objects are supposed to 
be inclined at small angles with respect to the line of sight (see 
Urry & Padovani, 1995 for a review), so the broad-line profiles 
will be significantly narrowed if disk-like BLRs are present per- 
pendicular to the jets, which can be used to test the geometry 
of BLRs in BL Lac objects. In this paper, we use the observed 
broad emission line widths and black hole masses derived from 
host galaxy luminosity to test the geometry of BLRs in BL Lac 
objects if the black hole masses can be estimated by indepen- 
dent method. The cosmological parameters £Im = 0.3, = 0.7, 
and Ho = 70 km s -1 Mpc" 1 have been adopted in this paper. 

2. BLACK HOLE MASS 

In order to estimate the central black hole masses of BL Lac 
objects, we use the relation between host galaxy absolute mag- 
nitude Mr at /?-band and black hole mass Mbh proposed by 
McLure & Dunlop (2004) 

log lo (M bh /M ) = -0.50(±0.02)Mr-2.75(±0.53). (1) 
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There are different surveys on host galaxies of BL Lac objects 
(e.g., Urry et al. (2000); Pursimo et al. (2002); Nilsson et 
al. (2003)). In this paper, we search the literature for all BL 
Lac objects with both measured host galaxies and broad emis- 
sion line profiles. As only a small fraction of BL Lac objects 
exhibit broad-line emission, we finally obtain a sample of 16 
BL Lac objects. All data collected are list in Table 1, and 
the derived parameters of these sources are listed in Table 2. 
The apparent magnitudes of the host galaxies at /?-band listed 
in Table 1 are galactic extinction and TT-corrected. Only up- 
per limits on host galaxy luminosity are available for 9 sources 
in this sample. We estimate the central black hole masses of 
these BL Lac objects from their host galaxy luminosity, and 
their broad-line profiles are used to explore their BLR geome- 
try. The black hole masses of two sources in our sample have 
been measured from the stellar dispersion velocity a, which 
gives log(M bh /M )=8.65(O521-365), 8.51(1807+698) (Barth, 
Ho, & Sargent 2003), and 8.90(1807+698) (Falomo, Koti- 
lainen, & Treves 2002), respectively. The black hole masses of 
these two sources derived from its host galaxy luminosity are 
log(M bh /M Q )=8.56(0521-365), 8.88(1807+698) (see Table 2). 
We can conclude that the black hole masses derived from host 
galaxy luminosity are quite reliable. 

3. BLR SIZES OF BL LAC OBJECTS 

The BLR size can be derived from the full-width half- 
maximum (FWHM) of the broad line 

GM hh 



Rblr = 



/ 2 ^FWHM 



(2) 



if the central black hole mass M bh is available. The correction 
factor / = 1/(2 sin/) for a pure disk-like BLR of which the axis 
is inclined to the line of sight at an angle i (McLure & Dunlop 
2001), while / = V3/2 for the clouds moving at random incli- 
nations (Wandel, Peterson, & Malkan 1999; Kaspi et al. 2000). 
Eleven sources in our sample have measured FWHM for broad- 
line Mg II, and four for Ha and another for H/3, respectively. 
The broad-line Mg II is expected to be produced in the same re- 
gion of the BLR, so we take VFWHivi(Mg II)=Vfwhm(H/3) for the 
source of which only the line width of H/3 is available. Using 
the black hole mass Mbh derived from the host galaxy luminos- 
ity and velocity Vfwhm, we can calculate the BLR size by using 
Eq. (2) on the assumption of isotropic motion of clouds in the 
BLR, i.e., / = vj/2 is adopted. We convert the size Rblr(Hcx) 
to /?blr(H/3) using relation (2) 

/WHa) = 1.19(±0.23)/?blr(H/3) + 13(±19) lt-day (3) 

in Kaspi et al. (2000) for the sources of which only the line 
widths of Ha are available. 

4. THE IONIZING LUMINOSITY 

BL Lac objects exhibit non-thermal continuum emission, 
which is believed to be dominated by the emission from the 
jets moving at relativistic speed at small angles with respect to 
the line of sight (Blandford & Rees 1978). The observed opti- 
cal continuum emission is a mixture of the emission from the 
disk and beamed emission from the jet, which prevents us from 
measuring ionizing luminosity directly from observed optical 
continuum emission for BL Lac objects (see Urry & Padovani, 
1995 for a review and references therein). Narrow-line emis- 
sion is suggested to be a good tracer for ionizing luminosity of 



radio-loud AGNs (Rawlings & Saunders 1991). We use the lu- 
minosity of the narrow-line [O II] at 3727 A to estimate the op- 
tical ionizing luminosity. For those the emission data of [O II] 
being unavailable, we estimate the [O II] emission from other 
emission lines adopting the line ratio proposed by Francis et al. 
(1991). The measured equivalent widths of narrow-line [Oil] 
corrected to the source frame corresponding to the measured 
continuum emission are usually around a few A or less than 1 
A for the BL Lac objects of which the narrow-line emission is 
detected (e.g., Stickel, Fried, & Kiihr (1989; 1993); Lawrence 
et al. (1996); also see Table 1 for the sources in our sample). As 
the observed optical continuum emission may be dominated by 
the beamed jet emission, the equivalent width EW; on of narrow- 
line [O II] corresponding to the ionizing continuum emission 
from the accretion disk should be larger than EW b s measured 
directly from the observed spectrum. The narrow-line emis- 
sion is isotropic and independent of the jet orientation, and we 
can therefore use the narrow-line luminosity L[oii] to estimate 
the ionizing luminosity at 3000 A for these BL Lac objects. 
Here we use the equivalent width of the line [O II]: EW; on =10 
A, corresponding to the ionizing continuum emission, which is 
suggest by Willott et al. (1999) to estimate ionizing luminosity 
of radio-loud AGNs. 

Using the same sample and the BLR sizes derived by Kaspi 
et al. (2000), McLure & Jarvis (2002) obtained a correlation 
between /? b ™r and the monochromatic continuum luminosity at 
3000 A 

/ XT \ 0.47+0.05 

« = (25.2 ± 3.0)(^) , (4) 

which is useful especially for the sources at high redshifts. In 
order to avoid confusion with the BLR sizes derived from Eq. 
(2), we use /? B ™ R to represent the BLR sizes derived from the 
empirical relation (4). 

5. RESULTS AND DISCUSSION 

We plot the relation between the ionizing luminosity 
AL A (3000) at 3000 A and BLR size R BL r in Fig. 1. The BLR 
size R-blr is derived from the width of broad line Mg II and 
black hole mass Mbh on the assumption of isotropic motion of 
clouds in the BLR, i.e., / = V3/2 is adopted in Eq. (2). We 
can also derive the BLR size R b "lr using the empirical relation 
(4). If the motion of BLR clouds is indeed isotropic, one may 
expect similar BLR sizes derived by these two different meth- 
ods. However, it is found that the sizes of BLRs Rblr in all 
these sources are ^2-3 orders of magnitude larger than R b "lr 
expected by relation (4). 

We note that only upper limits on the masses of the black 
holes in 9 of all 16 sources. The BLR sizes may be over- 
estimated for these 9 sources, because of the BLR sizes Rblr 
being derived from the line widths and black hole masses 
(see Eq. 2). The black hole masses estimated from the host 
galaxy luminosity for these 16 sources are in the range of 
~ 10 8 6 " 10 4 M Q (see Table 1). The deviations of the BLR sizes 
Rblr from /? B ™j> expected by relation (4) cannot be solely at- 
tributed to the overestimate of black hole masses for those nine 
sources with upper limits on galaxy luminosity, unless the black 
hole masses have been overestimated by 2-3 orders of magni- 
tude, i.e., the realistic black hole masses should be in the range 
of - 10 6 " 8 Mq for these sources, which seems impossible. It 
will be more difficult to attribute such deviations to the overes- 
timate of black hole masses for those 7 sources with well mea- 
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sured host galaxy luminosity. The black hole mass of the source 
1807 + 698 has been measured from its stellar dispersion veloc- 
ity a (Falomo, Kotilainen, & Treves 2002; Barth, Ho, & Sar- 
gent 2003), which is consistent with our estimate of the black 
hole mass 10 8 ,88 M . For this source, its BLR size Rblr derived 
from relation (2) is about three orders of magnitude higher than 
^blr predicted by relation (4) between R\™i\ and ALa(3000). 
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FIG. 1 . — The relation between the BLR size Rblr and the ionizing luminos- 
ity AL A (3000) at 3000 A derived from narrow-line luminosity L[o n l- Tne ^ ne 
represents the correlation Rg"^ - AL A (3000) defined by Seyfert 1 galaxies and 
quasars (McLure & Jarvis 2002). The source names are labelled in the plot. 

The ionizing luminosity at 3000 A is derived from the lumi- 
nosity of narrow-line [O II] adopting EW; on =10 A (Willott et al. 
1999), which is in general consistent with the optical spectro- 
scopic observations on the radio quasars selected from the Mo- 
longlo Quasar Sample (MQS), of which the equivalent widths 
of [O II] in the source frame are in the range from less than 1 
A to more than 100 A with an average of 14.7 A (Baker et al. 
1999). However, using a single value EW; on =10 A is still a 
rough estimate, and may induce uncertainties on the estimates 
of ionizing luminosity. In Fig. 1, we find that the deviations of 
Rblr from R b "lr w ^ not ^ e s °l ve d on ine isotropic BLR geome- 
try assumption even if a rather small EW; on =l A is adopted, i.e., 
the ionizing luminosity is an order of magnitude higher than the 
present values. In Table 1, we have listed the measured equiv- 
alent widths of [O II] in the source frame for the sources in our 
present sample, which are in 0.2-9.2 A. Considering the op- 
tical continuum emission contributed by the beamed emission 
from the jets in these BL Lac objects, their equivalent widths 
(EW; on ) corresponding to the ionizing optical continuum emis- 
sion should be larger than the measured values listed in Table 1 . 
This implies that the uncertainties on the estimates of the ioniz- 
ing optical continuum luminosity would not change the conclu- 
sion on this point. If photo-ionization of the gases in narrow- 
line region by radiative shocks driven by the radio source is 
important (e.g., Inskip et al. 2002), the central ionizing lumi- 



nosity should be lower than the present values, which would 
lead to even larger deviations. 

Wang & Zhang (2003) found that the BLR sizes of dwarf 
AGNs are systematically larger than the prediction of Rblr - 
XL\ correlation defined by Seyfert galaxies and quasars (Kaspi 
et al. 2000). They suggested that the flat ionizing spectra are in 
these dwarf AGNs as predicted by advection dominated accre- 
tion flow (ADAF) models, and the BLRs in those dwarf AGNs 
have lower ionization or/and lower density than those in Seyfert 
1 galaxies and quasars of Kaspi et al. (2000)'s sample. The 
sources in our present sample have brighter ionizing luminos- 
ity than those dwarf AGNs. As most BL Lac objects do not 
exhibit any line emission, all these 16 BL Lac objects have 
measured broad-line profiles and narrow-line emission have rel- 
atively high ionizing luminosity amongst all BL Lac objects 
(Cao 2002). There is a critical accretion rate tn crll , and an ADAF 
can exist only if accreting at a rate m < m cr ; t (e.g., Mahadevan, 
1997), which leads to an upper limit on optical continuum emis- 
sion from an ADAF for a given black hole mass (Cao 2002). 
Cao (2003) calculated optical spectra of ADAF+SD(standard 
disk) systems and compared them with the observed spectra, 
which suggests that only ADAFs themselves are unable to pro- 
duce such bright ionizing optical continuum emission and stan- 
dard thin disks should be present at least in the outer regions 
of the disks for these BL Lac objects. The ionizing photos are 
therefore mainly from the standard thin disk regions in these 
BL Lac objects, unlike the dwarf AGNs considered by Wang 
& Zhang (2003). The luminosity ALa(3000) of the sources in 
the sample used to derive the correlation Rblr - ALa (3000) is in 
the range of - 10 34 " 39 W (McLure & Jarvis 2002). The ionizing 
luminosity of the BL Lac objects in present sample at 3000 A 
are in the range of ~ 10 35 8-38 W, which is in the similar range 
as their sample. The Rblr- ALa (3000) relation derived from 
a sample of quasars and Seyfert galaxies should be valid for 
these BL Lac objects, unless the physics of BLRs in these BL 
Lac objects is significantly different from that in the sources of 
the sample considered by McLure & Jarvis (2002). 

We derive the BLR sizes of these BL Lac objects assuming 
the motion of clouds to be isotropic, which may not be the case 
in these sources. Our estimate of BLR sizes may strongly be 
overestimated due to anisotropic cloud motion, if the velocity 
component projected to the line of sight is only a small fraction 
of its real velocity. A most likely candidate for such anisotropic 
motion of clouds is the clouds orbiting in a disk-like BLR (e.g., 
McLure & Dunlop, 2001). For the clouds orbiting in a disk-like 
BLR, the correction factor / in Eq. (2) is 1 /(2sin0, where i is 
the angle of the axis inclined to the line of sight (McLure & 
Dunlop 2001). If this is the case, we can estimate the inclina- 
tion angle i of these BL Lac objects assuming they indeed to 
obey the correlation Rblr- ALa(3000) suggested by McLure 
& Jarvis (2002), i.e., the value of / is estimated by letting 
Rblr = ^blr- The derived results are listed in Table 1. We find 
that the inclination angles are around <~ 2°- 12° for these BL 
Lac objects. There is evidence that the velocity field of BLR is 
better described by a combination of a random isotropic com- 
ponent, with characteristic velocity V r , and a component only 
in the plane of the disk, with characteristic velocity V p (Wills & 
Browne 1986). In this case, the observed FWHM will be given 
by 

VFWHM = 2(V r 2 + V p 2 sin 2 1/2 (5) 

(McLure & Dunlop 2001), so / = 0.5[(y r /y p ) 2 + sin 2 /]- 1 / 2 . If 
the random isotropic component is important, i.e., V r is compa- 
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rable with V p , then the term (V r /V ? ) 2 cannot be neglected and 
the derived inclined angle of the disk axis will be less than that 
listed in Table 2. This is in general consistent with the unifi- 
cation schemes that the jets of BL Lac objects are inclined at 
small angles to the line of sight. 
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Table 1 



Data of BL Lac objects 



Source 


Redshift 


mff(host) 


References 


FWHM(Mg 11) 


References 


log \L X 


EW of [O n] 


References 


0138-097 


0.733 


>18.38 


1 


4842 


4 


31.23, 


0.52 


4 


0521-365 


0.055 


14.35 


1 


3000"" 


5 


36.30" 




11 


0715-259 


0.465 


16.76 


1 


5200 


6 


36.70 


1.3 


6 


0820+225 


0.951 


>19.38 


1 


2995 


7 


35.84 e 




7 


0823+033 


0.506 


> 19.04 


1 


5455 


7 


36.53 


0.2 


7 


0851+202 


0.306 


18.44 


2 


2635 c 


8 


37.28 f 




8 


0954+658 


0.367 


>18.85 


1 


2079" 


9 


36.72 


0.22 


4 


1144-379 


1.048 


> 19.93 


1 


2492 


8 


37.13° 




8 


1308+326 


0.996 


>18.36 


2,3 


4016 


7 


37.45° 




7 


1538+149 


0.605 


18.73 


1 


2411 


7 


37.15 


0.44 


7 


1803+784 


0.684 


>19.15 


1 


3082 


4 


37.41 r 




9 


1807+698 


0.051 


13.54 


1 


1326" 


3 


36.50 


2.5 


7 


1823+568 


0.664 


18.57 


1 


3952 


4 


37.54 


1.9 


7 


2131-021 


1.285 


> 18.50 


1 


3602 


4 


37.73 


9.2 


4 


2200+420 


0.069 


14.55 


1 


4260" 


10 


35.46 f 




7 


2240-260 


0.774 


>20.22 


1 


2753 


7 


38.06 


1.4 


7 



"FWHM of Ha. 

b The profile of broad-line Ha in this source is asymmetric with a red wing of FWHM ~ 3000 km s -1 and a blue wing of FWHM ~ 
1500 km s -1 (Scarpa, Falomo, & Pian 1995). Here, we conservatively take FWHM=3000 km s -1 . 

c FWHMofH/3, 

[l the equivalent width of [O II] is not given, 

"the flux of [O II] is converted from Mg II by using the ratio given in Francis et al. (1991), 
r the flux of [O It] is converted from [O 111] by using the ratio given in Francis et al. (1991), 



Note. — Column (1): source name; Column (2): redshift; Column (3): galactic extinction and AT-corrected fi-band magnitude of the host 
galaxy; Column (4): references for R-band magnitude of the host galaxy; Column (5): FWHM of broad-line Mgll (km s~'); Column (6): 
references for FWHM(Mgll); Column (7): ionizing luminosity at 3000 (W) estimated from Ljoirj; Column (8): measured equivalent width 
of [O II] corrected to the source frame (); Column (9): references for Z4011] ■ 



References. — (1) Urry et al. (2000); (2) Pursimo et al. (2002); (3) Nilsson et al. (2003); (4) Rector & Stocke (2001); (5) Scarpa, Falomo, 
& Pian (1995); (6) Carangelo et al. (2003); (7) Stickel, Fried, & Kiihr (1993); (8) Stickel, Fried, & Ruhr (1989); (9) Lawrence et al. (1996); 
(10) Corbett et al. (1996); (11) Tadhunter et al. (1993). 
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Table 2 



Derived parameters of BL Lac objects 



Source 


log M bh /M Q 


log Rblr 


/ 


i 


0138-097 


<9.70 


<3.16 


<5.8 


>4.9° 


0521-365 


8.56 


2.34 


4.2 


6.9° 


0715-259 


9.91 


3.31 


9.2 


3.1° 


0820+225 


<9.55 


<3.43 


<16.8 


>1.7° 


0823+033 


<8.88 


<2.24 


<2.9 


>9.8° 


0851+202 


<8.78 


<2.78 


<3.6 


>7.9° 


0954+658 


<8.56 


<2.68 


<4.8 


>5.9° 


1144-379 


<9.40 


<3.44 


<8.5 


>3.4° 


1308+326 


<10.12 


<3.75 


<10.1 


>2.8° 


1538+149 


9.27 


3.34 


7.4 


3.8° 


1803+784 


<9.22 


<3.08 


<4.8 


>6.0° 


1807+698 


8.88 


3.46 


12.2 


2.3° 


1823+568 


9.47 


3.11 


4.7 


6.2° 


2131-021 


<10.39 


<4.11 


<13.3 


>2.2° 


2200+420 


8.71 


2.18 


5.5 


5.2° 


2240-260 


<8.85 


<2.81 


<2.5 


>11.7° 



Note. — Column (1): source name; Column (2): black hole mass; 
Column (3): BLR size for Mg II (light-day); Column (4): derived cor- 
rection factor / for BLR geometry; Column (5): inclination angle of 
the jet with respect to the line of sight. 



